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; MTHSATAN is a function watch returns the g Testing point arc tangent 
; value (in radians) of its s Inge precision floating point orgpent 
; MTHSATAN2 is two arguments floating point arctangent. The call i 
; standard call-by-reference. 
; MTHSATAN_R4 is a special routine which is the same as MTHSATAN 
; except a faster non-standard JSB call is used with the argument in 
; ; RO and no registers are saved. 


:  MTHSATAND is a function which returns the floating point arc tangent 
value (in degrees) of its single precision floating point argument. 
: NTHSATAND D2 is two arguments floating point arctangent. The call is 
; standard call-by-reference. 
: MTHSATAND_R4 is a special routine which is the same as MTHSATAND 
; except a faster non-standard JSB call is used with the argument in 
; RO and no registers are saved. 
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Add "_"' to the PSECT directive. JBS 21-DEC-7 

Define all global symbols using .EXTRN. JBS 19-JUN-1979 

Added degree entry points. RNH 15=MAR-1981 
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16-SEP-1984 01:04:23 YAX/VMS Macro V04-00 
=SEP=1984 71:95:87 MTHRTL.SRCIJMTHATAN.MAR; 1 


-SBTTL MTHSATAN = Standard Single Precision Floating Arc Tangent 


Precision Flo 


++ 

; ; FUNCTIONAL DESCRIPTION: 

; : ATAN - single precision floating point function 
: : ATAN is computed using the following steps: 


1. If X > 5 then 
a. tet u = 1/X. 
b. Compute ATAN(W) = = WeP(Wee2), where P is a polynomial of 


C. gett ATAN(X) = s0'le 00 - ATAN(W) 


me 3F 
a. Obtain: xi “by table. Look-up. 
b. Compute X = MHI)/(1 + X*KXHI). 


XHI 
C. compute ArANCD) = 1*P(7**2), where P is a polynomial of 


2. 
d. oe rese ATAN(XHI) by table look-up. ATAN(XHI) will have 
two parts - the high order bits, ATAN_XHI_HI, and the low 


order bits, ATAN_X I 
2 ATEN _XHI HI + (ATAN_XHI_LO + ATAN(Z)). 


e. ree ATAN(X) = 
3. If O = < 3/32 then 
a. lonmeke “ATANOXD. = X + X#Q(X*"2), where Q is a polynomial 
of egeree 2. 
4. 1f x < Ot 
* Co meme nv = te Xi) using steps 1 to 3. 


Set ATAN(X) = -¥ 


CALLING SEQUENCE: 
Arctangent.wf.v = MTHSATAN(x.rf.r) 


; INPUT PARAMETERS: 


; define longword multiplier 
; x is an angle in radians 


LONG = 4 

x = 1 * LONG 
IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 

VALUE: floating arctangent angle of the argument 
IMPLICIT OUTPUTS: none 
SIDE EFFECTS: 
Signals: none 
NOTE: This procedure disables floating point underflow and integer 


overflow, causes no floating overflow or other arithmetic traps, and 
preserves enables across the call. 


ns 
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2-005 MTHSATAN = Standard Single Precision Flo 6=SEP-1984 11:20:51 CMTHRTL.SRCJMTHATAN.MAR; 1 (6) | s 
1BC 399 ;==- A 
ee 88 
401c O1BC 4 : -ENTRY MTHSATAN, ACMASK ; standard call-by-reference entry . 
1BE 40 i; disable DV (and FU), enable IV 
18 404 MTHSFLAG_JACKET ; flag that this is a jacket procedure r 
6D 00000000'GF 9E 018 MOVAB G*MTHSSJACKET_HND, (FP) A 
1C ; set handler address to jacket 
i ; handler A 
1¢ 405 ; in case of an error in special JSB A 
1¢5 4 $ ; routine A 
50 04 8¢ 50 01¢5 40 MOVF ax(AP), RO ; RO = arg ; A 
2 10 01C9 408 BSBB = MTHSATAN_R4 ; call special ATAN rountine A 
04 pice 409 RET 3 return = result in RO I 
1cC = 410 : 
L 
L 
L 
L 
L 
L 
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~-SBTTL MTHSATAN2 = Standard floating Arctangent With 2 Arguments 


39 
; FUNCTIONAL DESCRIPTION: 


If Y > 
If Y< 


CALLING SEQUENCE: 


ATAN2(X,Y) 
(X,Y 
5, ATAN2(X,Y 


AO ODONA UE WIN $$ O ODNOU EW" OOONAUS UMN oF 


> ¥ <0 and X < 0 and X/Y 


3% 91:06:63 AX/VMS Macro V04-00 


MTHRTL.SRCIMTHATAN.MAR; 1 


; ATAN2 = single precision floating point function 
; ATAN2(X,Y) is computed as following: 
If Y= or X/Y¥ > 2##25 


and X/Y =< 
and X/Y =< 


~~ 


Arctangent2.wf.v = MTHSATAN2(x.rf.r, y.rf.r) 


; INPUT PARAMETERS: 


x = 1 * LONG 
y = 2 * LONG 


ENTRY MTHSATAN2 ,ACMASK 
MTHSFLAG_JACKET 


MOVAB G*MTHSSJACKET_HND, (FP) 


MOVF ax(AP), RO 
MOVF ay(AP), R1 


5 3 
; Test if Y = 0 or X/Y > 2*#25 


INF 
BICW #*x807F, RO, Ré 
ft #*xX807F, R1, R 


. Re 
R2, #26*128 
BGTR _INE 


: Test if Y>Oory<0 


R1 
BGTR . PLUS 
BGEQ A1PLUS 


=< 22225 


; x is the first argument 
; y is the second argument 


; SIDE EFFECTS: See description of MTHSATAN 


standard call-by-reference entry 


disable DV (and FU), enable IV 


flag that this is jacket procedure 


set handler address to jacket 
handler 


; in case of an error in special JSB routine 
arg 


; RO = 
; R1 = arg2 
branch to INF if Y = 0 
R2 = exponent (X) 
R35 = exponent (Y) 
R2 = exponent(X) = exponent (yY) 
compare RS with 26 
if X/Y > 2**25, branch to INF 


test the sign of Y 
branch to AdPLUS if Y > 0 
test the sion of Xx 
branch to AIPLUS if X >= 0 


vz 
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2-005 ATHSATAND == Standard flost ing Arctangent grSEb a1 984 9:90:87 MTHRTL.SRCIJMTHATAN.MAR; 1 (3) 
B 1 1 464 BSBB MTHSATAN_R4D ; RO = ATAN(X/Y) 
50 OFDB4149 H 1 1 465 SUBF #LF_PI, RO 2 Re = =P] + ATAN(X/Y) 
04 : § RET 3; return 


6 

6 

‘ ; Y < 0 and X > O and X/Y =< 2825 
A1PLUS: 

7 

Ff 

7 


0 
Ht 18 1 BSBB MTHSATAN_R4D s Re = ATAN(X/Y) 
50 OFDB4149 8F 4 ¢ ADDF LF PI, RO ; RO = PI + ATAN(X/Y) 
04 RET 3 return 
> Y¥ > 0 and X/Y¥ =< 20825 
A2PLUS: 
27 10 BSBB MTHSATAN_R4D ; RO = ATAN(X/Y) 
04 RET ; return 


OOOCCCOCOCOCOCOCOOOCOOOCOO 


Sel tek delta dee del ee ee ee ee ae ae ee et at et 


74 
g? 
§ 
78 
3 
H : Y = 0 or KX/Y¥ > 22225 
0 § : 
0 83 INF: 
50 B5 0 84 TSTW RO ; test the sign of X 
OA 14 0 85 BGTR 1$ : branch if x > 0 
— 2 © 6 BEQL 2$ : branch if K = 0 
50 OFDBCOC9 BF 50 0 ¥ MOVF § #LF_MPI_OVER_2, RO : RO = ATAN(X/Y) = =PI/2 
04 8 Hy RET 3 return 
50 OFDB4OC9 BF 50 0 90 1$: MOVF § #LF_PI_OVER_2, RO ; RO = ATAN(X/Y) = PI/2 
04 0 491 RET 3 return 
0 $38 
0 493 5+ ; 
0 494 ; Here if X = 0 and Y = 0. Signal INVALID ARG TO MATH LIBRARY 
02 495 ; as a SEVERE error. 
02 496 ;- 
Bs 497 
50 01 OF 79 498 2$: ASHQ #15, #1, RO ; RO/R1 = reserved operand which 
022 499 ; is copied to CHFSL_MCH_SAVRO/R1 
022 500 3; so a handler can fixup if wants 
0 501 3; to continue 
7E 00° 8F 9A 0 288 MOVZBL #MTHSK_INVARGMAT, -(SP) ; code for INVALID ARG TO MATH LIBRARY 
00000000'GF 61 FB Bs 50 CALLS #1, G*ATHSSSIGNAL 3 signal SEVERE error 
04 0 504 RET : return if handler continues 
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2-005 ATHSATAN. RG = Special ATAN routine b-SE b=} 984 9:90:87 MTHRTL.SRCIJMTHATAN.MAR; 1 = (8) | 
rn ~SBTTL MTHSATAN_R4 = Special ATAN routine | 
rn : Special ATAN = used by the standard routine, and directly. 
A 910 ; CALLING SEQUENCES: 
A 12 save anything needed in RO:R4 
A 1 ; MOV ase ; input in RO 
A 513; JSB MTHSATAN_R4 
2 13 : return with result in RO 
A 1 ; Note: This routine is written to avoid causing ony integer overflows, 
A 7; floating overflows, or floating underflows or divide by 0 conditions, 
8 . 18 3 whether enabled or not. 
A 0 ; REGISTERS USED: 
A 521; RO - ¥ Lost ing argument then result 
A 2 ¢ ; RO:R3 = POLY 
023A 3 R4 = Pointer into ATAN_TABLE 
0 - 2 4 
023A 5 
023A 5 § 
023A 527 MTHSATAN_R4D 
50 51 46 O23A 528 V R1, RO : for our own use only! 
8 D 529 MTHSATAN_R4:: : Special ATAN routine 
50 53 023D 530 | TSTF RO ; R4 = X = argument ; 
7% «©6119 Q ef : 1 BLSS NEG_ARG ; Branch to negative argument Logic 
p 41 338 ; Argument is positive 
54 50 3€CO BF Ad 0241 535 ° SUBW3 #*X3ECO, RO, R4 ; Argument is less than 3/32 
45 19 0247 536 BLSS SMAL ; branch to small argument logic 
54 O36F 8F Bi 0249 537 CMPW #*X0356F, R4 ; Argument is greater that 11, 
41 19 8 3 238 BLSS LARGE _ARG ; branch to Large argument logic 
0 23 379 : Logic for positive medium sized arguments. Get pointer into ATAN_TABLE. 
54 54 FC BF 9 6 2 235 : ROTL #-4, R4, R4 ; R4 = index into MTHS$AB_ATAN table 
54 FFFFFFOO 8F CA 8 5 4 BICL #-256, ; zero high order bits of index 
54 FDA3 CF44 90 5 544 MOVB MTHS$S$AB_ATANC(R4], R4 ; R4 = offset into ATAN_TABLE 
54 FDDS CF44 DE 0 ° 22 MOVAL ATAN_TABLECR4], R4 : R4 = pointer to XHI 
6 i : Compute Z 
51 84 09 6 Pr , MOVL. (R4)+, R1 : RI = XHI 
52 50 1 4 68 50 MULF3 R11, RO, R2 : Re = X*XHI 
52 8 3 f 551 ADDF #1, R 3 R2 = 1 + X*XHI 
50 1 4 g 236 SUBF R1, R ; RO = ¥ = XHI 
50 2 46 4 37 DIVF R2, RO 3; RO = Z = (X = XHID/(1 + X*XHI) 
i 33 > Evaluate Z*P(Z*#2) 
7 50 00 7 25 . MOVL RO, -(SP) : Push Z onto the stack 
5p 33 4&4 78 8 MULF RO, R > RO = Z#*2 
FE60 CF 36002 Ss 50s ze 39 POLYF RO, WATANLENT~1, ATANTABT pinkenh 
50 Hi GG be 2) MULF (SP)+, RO ; nd : ATAN(Z) = Z*P(Z7**#2) 
50 84) 40 7 562 ADDF = (R4) +, RO ; RO = ATAN_XHI_LO + ATAN(Z) 


—z 
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MTHSATAN ; Floating Point Arc Tangent Functions 16-SEP-1984 01:04:23 VAX/VMS Macro v04-00 Pa 14 
2-005 ATHSATAN.RG = Special ATAN rout ine BoSEP-19R6 $1:56587 FATTY cacdortintageman:1 29° 1d, 
50 64 40 A 63 ADDF (R4), RO ; RO = ATAN(X) = ATAN_XHI_HI + 
D 564 : (ATAN_XHI_LO +~ATAN(Z)) 
05 4 68 RSB 3; Return 
E 6 
0096 = 31 ‘ oe SMALL: BRW SMALL_ARG ; Dummy Label used to avoid adding 
8 : an extra insrtuction in the 
4) oy ; medium argument logic 
4) ie ; Large positive argument logic. 
91 574° 
91 575 LARGE_ARG: 
eee oak My eae St MULES Rao Re RG’ izces 
Fe41 CF «60025055 ons Hi POLYF RO, #ATANLENT-1, ATANTABI coptiae 
50 54 44 6 A 80 MULF R4, RO : RO = ATAN(W) = -WeP(Wee2) 
50 BD2EB438 BF 40 02A6 581 ADDF  #Lf_PI_OVER_2.LO, RO : 
50 OFDB4OC9 8F 40 O2AD 8¢ ADDF #LFIPI“OVERT2-HI, RO =: RO = ATAN(X) = PI/2 = ATAN(W) 
05 8 B4 28 RSB 3; Return 
0 Be A ; 
0285 286 : Logic for negative arguments 
oops ea | 
6 B5 589 NEG_ARG: 
54 50 BECO 8F A3 0285 590 SUBW3 #*XBECO, RO, R4 ; Argument is less than 3/32 
A 19 Q2BB 591 BLSS SmALL ARG ; branch to small argument logic 
54 O36F 8F 81 O2BD 238 CMPW #*XO35F, R4 ; Argument is greater than 11, 
7 WW 8 Cf 27 BLSS N_LARGE_ARG ; branch to large argument logic 
8 $2 232 : Logic for negative medium sized arguments. Get index into ATAN_TABLE. 
54 54 FC or 9C 8 C4 299 . ROTL #-4, R4, RS ; R4 = index into MTH$$AB_ATAN table 
54 FFFFFFOO 8F CA C9 «86598 BICL -256, R ; clear high order (unused) bits of ind 
54 FD2F CF44 90 8 dO 599 MOVB MTHS$SAB_ATANCR4], R4 ; R4 = offset into ATAN_TABLE 
54 FD61 CF44_~ «ODE 43 $99 MOVAL ATAN_TABLECR4], R4 ; R4 = pointer to XHI 
BEE Gig | commute 2 
51 8 DO O2DC 604 ° MOVL. (R4)+, R1 : R1 = XHI 
52 5051 43 DF 605 MULF3 R1, RO, R2 > R2 = X*XHI 
5 98 ¢ 4 + O06 SUBF R2, #1, R 3 R2 = 1 = X*XHI] = 1 + Xe(=KXHI) 
40 E 6 ADDF ° 3; RO = X + XHI = X = (-XHI) 
50 2 46 eA on8 : DIVF R2, RO : RO = 
ED 610 : Evaluate 2*P(2*#2) 
7E 50 00 OED 61¢ ; MOVL RO, -(SP) : Push Z onto the stack 
50 50 44 F 61 MULF . >; RO = Z**2 
FOEB CF §=6(02Ss«50Ss«*S55 3 614 POLYF RO, #ATANLENI-1, ATANTABI soi Seaaiaihd 
50 gE 44 02F9 616 MULF (SP)+, RO : RO = ATAN(Z) = 2*P(Z##2) 
4 $s FC 61 SUBF (R4)+, RO ; RQ = ATAN_XHI_LO + ATAN(Z) 
50 64 4 FF o18 SUBF (R4), RO : RO = ATANTX) = ATAN_XHI_ HI + 
02 61 ; (ATAN_XHI_LO +7ATAR(Z)) 


uz 
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2-005 ATHSATANTRG © Special ATAN routine 6-S Ep=19 1987 9 £96: $7 MTHRTL.SR ar MTHATAN.MAR; 1 . 
05 9 RSB ; Return 
§ ; Logic for Large negative arguments 
i e 
é 5 N LARGE _ARG: 
54 0000C080 ef 0 47 § “DIVF RO, #-1, R4 3 RG = W = 1/IX! 
50 4 4 $2 B 6 MULF R4, RO > RO = Wee2 
Focr CF 4602 «660~Co 5 F 628 POLYF RO, #ATANLEN1=1, ATANTABI 
13 629 7: RO = P(Wee2) 
50 4 44 1 630 MULF R4, RO 3 RO = ATAN(W) = WeP(Wee2) 
50 BD2EB4 3B F 6 18 6 1 SUBF § #LF_PI_OVER_2_L0, R : 
OFDB40C9 BF 4 1F ; SUBF § #LF-PI“OVERTO“HI. R : RO = ATAN(X) = ATAN(W) = PI/2 
0 § ° ? RSB 3; Return 
a5 5 OE smatt t Logi 
3 Small argument logic. 
327 $39 : 
; 7 o98 
7 639 SMALL_ARG: 
54 DO 0327 640 MOVL RO, 3 R4 = prquaent = X 
50 8000 8F AA O32A 641 BICW FK8000, RO ; RO = 
50 3a00 8F 861 O32F S66 CMPW #*x3A00, RO ; Compare "24-13 to {X! 
19 0334 64 BLSS 1$ 3 Branch to polynomial evaluation 
50 54 D0 0336 644 MOVL R4, RO 3 No POLY needed. Answer = X 
05 0339 645 RSB 3 
033A 646 
50 50 44 QO33A 647 18: MULF RO, RO > RO = X**2 
FDAD CF 4602) «C«50s«S5s«O33D = 648 POLYF RO, #ATANLEN2=1, ATANTAB2 
034 649 3; RO = Q(X*##2) 
50 54 44 034 650 MULF R4, RO : RO = X*Q(X##2) 
50 54 40 0346 651 ADDF R4, RO ; RO = ATAN(X) = X + X#Q(X##2) 
05 0349 626 RSB ; Return 
034A 65 
034A 654 
034A 655 
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-SBTTL MTHSATAND = Standard Single Precision Floating Arc Tangent 


FUNCTIONAL DESCRIPTION: 
ATAND = single precision floating point function 
ATAND is computed using the following steps: 


1. If X ; > i, een 
a. Let 
b. tonpute ATAND(W) = WeP(We*2), where P is a polynomial of 


de 

3 ef at ATAND(X) = pie - ATAND(W) 
/32 =< X =< 11 ¢t 

Obtain XHI by tabte’ look-up. 

Compute Z = (X = XHI)/(1 + X*XHI). 

C. enone ATAND(Z) = Z*P(Z* 22). where P is a polynomial of 


d. oe rose ATAND(XHI) by table Look-up. ATAND(XHI) will have 
two parts - the hi at wal bits, ATAND_XHI_HI, and the low 
order bits, ATAND_ LO. 

e. Compute ATAND(X) = ATAND _XHI_HI + (ATAND_XHI_LO + ATAND(Z)). 

3. If 0 =< X < 3/32 then 

a. Compute ATAND XD = 64*X + X*#Q(X"#2), where Q is a polynomial 

of degree 2. 


2. 
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:* 
: 
; 
: 
; 
; 
: 
= 
: 
Py 
. 
Py 
. 
: 
. 
: 
; 
Py 
rd 
‘ 
Py 
7 
Py 
: 
: 
. 
Py 
: 
Py 
. 
: 
: 
Py 
a 


4A 
4A 
4A 
4A 
4A 
4A 
4A 
4A 
4A 
4A 
034A 
4A 
GA 
034A 
GA 
034A 
034A 
034A 68 4. If x <0 then 
034A $8 a. Compute Y = oa iy Xi) using steps 1 to 3. 
034A 6 b. Set ATAND(X) = -¥ 
034A 68 
034A 68 
034A 689 ; CALLING SEQUENCE: 
034A 69 
034A 691 ; Arctangent.wf.v = MTHSATAND(x.rf.r) 
034A 69 
034A 693 ; INPUT PARAMETERS: 
been 694 
00000004 034A 695 LONG = 4 ; define longword multiplier 
00000004 Baca 696 x = 1 * LONG ; x is an angle in radians 
034A 698 ; IMPLICIT INPUTS: none 
034A 344 3 
Been £80 ; OUTPUT PARAMETERS: 
Baen me 3 VALUE: floating arctangent angle of the argument (in degrees) 
3h 4 : 3 IMPLICIT OUTPUTS: none 
re 7 $ 3; SIDE EFFECTS: 
4A 707; 
rn a 3 Signals: none 
4A 710 ; NOTE: This procedure disables floating point underflow and integer 
4A 711; overflow, causes no floating overflow or other arithmetic traps, and 
r+ ne 3 preserves enables across the call. 
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.ENTRY MTHSATAND, ACMASK 
MTHSFLAG_JACKET 
G*MTHSSJACKET_HND, (FP) 
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rat AX/VMS Macro V04-00 Page 
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; standard call-by-reference entry 
; disable DV (and FU), enable IV 
; flag that this is a jacket procedure 


; set handler address to jacket 
handler 


; in case of an error in special JSB 


; routine 
; RO 


= arg 
; call special ATAND coutone 


; return - result in 


PONWAN 


wz 


MTHSATAN ; Float 
2-005 ATHSATA 


ao 


oin 
Standard floating Arctangen 6-SEP-1984 MTHRTL.SRCJMTHATAN.MAR; 1 
7 = ~SBTTL MTHSATAND2 = Standard floating Arctangent With 2 Arguments 
; FUNCTIONAL DESCRIPTION: 
: ATAND2 = single precision floating point function 

4 

5 

6 
i 


SA 7 
: 7 
a 
ee 
A 
sh if ; ATAND2(X,Y) is computed as following: 
SA 735: If Y = 0 or X/¥ > 20825 ATAND2(X, ¥) = 90 * (sign x) 
0 SA 786: If Y > 0 and X/Y =< "505, A AND G(X.) = ATAND(K/Y) 
O34 if ; If ¥ < 0 and X/Y =< 20025; ATAN ND2(X,Y) = 180 * (sign X) + ATAND(X/Y) 
A 739: 
oe2h oo ; CALLING SEQUENCE: 
tn rs : Arctangent2.wf.v = MTHSATAND2(x.rf.r, y.rf.r) 
O34 ad : INPUT PARAMETERS: 
00000004 O35A 746 x = 1 * LONG ; x is the first argument 
00000008 Bae re y = 2 * LONG ; y is the second argument 
O35A 749 ; SIDE EFFECTS: See description of MTHSATAND 
035A 750; 
035A 751 ;:-=- 
Bish 
401C Q35A 754 -ENTRY MTHSATAND2 ,ACMASK ; standard call-by-reference entry 
035C 755 ; disable DV (and FU), enable Iv 
baec 756 MTHSFLAG_JACKET ; flag that this is jacket procedure 
6D O00000000'GF 9E b326 MOVAB G*“MTHSSJACKET_HND, (FP) 
036 ; set handler address to jacket 
baee 3; handler 
0363 ree ; in case of an error in special JSB 
363 «= 758 3 routine 
50 04 8C 50 0363 759 MOVF ax(AP), RO ; RO = argl 
51 08 BC 50 re 760 MOVF @y(AP), R1 ; R1 = arg2 
0366 166 : Test if Y = 0 or X/Y > 20#25 
35 13 $368 764 BEQL 3; bra nch to INFD if Y = 0 
26 50 7F or AB $9 765 BICW Be > RO. - 3 if exponent (X) 
5 51 7F 8F . 76 BICW a*yx Rp ; RS = exponent(yY) 
ze 53 A 79s 76 SUBW Re. ; R2 = exponent it? - exponent (Y) 
OD00 8F 52 8B) 7C poe CMPW #26*128 3; com pare “ with 26 
| if HI it: BGTR Ree : if K/Y > 2*25, branch to INFD 
; 7 : Test if Y>Oory<0 
51 B5 ng ’ TSTw = R11 ; test the sign of Y 
18 14 8 774 BGTR A2PLUSD 3; branch to ; PLUSD if y>oO 
5 83 ¥: ore TSTW 3; test the s sign ¢ 
OA 1 M ors BGEQ AIPLUSD : branch to A PLUSD if xX >= 0 
B 778: Y <0 and X <0 and X/Y =< 2##25 
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MTHSATAN 
2-005 


2 
00004434 8F 


50 
H 
50 «00004434 8F 
27 
50 
OA 
10 
50 0000C384 8F 
50 00004384 8F 
50 01 OF 
7E _00°8F 
00000000'GF 01 


oul ourw—@ 
Fo FOWSM 


9A 


paring Point Arc Tangent banat tone 1 
ATANDZ = Standard floating Arctangen 
B 779; 
8 780 6SBB MTHSATAND_&4D 
D 781 SUBF #LF_180, RO 
2 ¢ RET 
95 783; 
9 784; Y < 0 and X > O and X/Y =< 2##25 
9 785 ; 
1 eee $ A1PLUSD: 
95 7 BSBB MTHSATAND_R4D 
97 788 ADDF #LF_180, RO 
9E 789 RET 
te OST |v > © ond A/4 0e 20025 
5 > an =< grt 
oF 35 : 
039F 793 A2PLUSD: 
O39F 794 BSBB ss MTHS$ATAND_R4D 
3A1 795 RET 
a5 $89 3 yw 0 or xiv > 20025 
: 2 or > 2re 
O3A 798 ; 
O3A 799 INFD: 
3A 800 TSTW RO 
03A4 = 801 BGTR 1$ 
03A6 ane BEQL 2$ 
O3A8 80 MOVF  § #LF_M90, RO 
Baar 804 RET 
380 805 
0380 806 1$: MOVF #LF_90, RO 
038 807 RET 
o3e8 808 
038 809 ;+ 
0388 10 ; Here if X = 0 and Y = 0. 
tet) 11 ; as a SEVERE error. 
0388 \§ = 
tot 81 
038 814 28: ASHQ #15, #1, RO 
038C 815 
03B8C 816 
038C 836 817 
038C 818 MOVZBL  bige INVARGMAT, 
03C0 819 CALLS . =” ATHS$S1GNAL 
03C7 ~=—820 


wSP=1986 1: 


-(SP) 


26:87 


; we De 


AX/VMS M 
MTHRTL.SR 


v04- 


bo 4 


Re = ATAND(X/Y) 
RO = -180 + ATAND(X/Y) 
return 


RO = ATAND(X/Y) 
RO = 180 + ATAND(X/Y) 
return 


RO = ATAND(X/Y) 
return 


test th 
branch 


P--—-® 


return 


ATAND(X/Y) 
turn 


90 


Signal INVALID ARG TO MATH LIBRARY 


ro 00 
RCIMTHATAN. MAR: 1 


RO/R1 = reserved operand which 


be co180ed to CHF$ 
2 Sooeet can fixup 
to tent tay 


_MCH_SAVRO/R1 


f wants 


Page 


code for INVALID ARG TO MATH LIBRARY 


signal SEVERE error 


return if handler continues 


19 


(10) | 


vs 


MTHSATAN 
2-005 


-v 


D 12 | 
oint Arc Tangent Functions 1975601386 91:08:93 AX/VMS Macro Vv04-00 Page 20 
- Special ATAND routine 6-SEP=-1984 11:20:51 MTHRTL.SRCIMTHATAN.MAR; 1 (11) 


-SBTTL MTHSATAND_R4 = Special ATAND routine 
Special ATAND = used by the standard routine, and directly. 


CALLING SEQUENCES: 
neve anything needed in RO:R4 


JSB MTHSATAND_R4 


=z: 
— 
x= 
An- 
>o 


| 

; input in RO | 

| 

return with result in RO | 
| 


Note: This routine is written to avoid causing any integer overflows, 
floating overflows, or floating underflows or divide by 0 conditions, 
whether enabled or not. | 


REGISTERS USED: 
RO - Floating argument then result 
RO:R3_ = POLY 
R4 = Pointer into ATAND_TABLE 


Sete Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge 


50 51 46 R1, RO ; for our own use only! 
MTHSATAND_R4:: : Special ATAND routine 
50 53 TSTF RO ; R4 = X = argument . 
6D 19 BLSS NEG_ARGD ; Branch to negative argument logic 


4 

8 

9 

0 

4 

8 

9 

; 

§ MTHSATAND_R4D: 
4 

5 

8 . 

8 : Argument is positive 
0 
1 
2 
3 
4 
5 
6 
7 
8 
9 
0 
4 


BPRS FFA NA NNAAWWNAW AANA WWW AINA AWWW NAAN oe 


SOOOOOCOCOOOCOOOOOCOCOOCOSOOOCCOOSOSCCOCCSCCOCOCCOCOOOSOCOCOCOCOOCOOOOOCOOCOCoOOoOO —a 


54 50 36CO 8F A3 5 SUBW3 = #*X3ECO, RO, R4 ; Argument is less than 3/32, | 
19 5 BLSS asta 4 ; branch to small argument logic 

54 O36F 8F 861 5 CMPW #*XO036F, R4 ; Argument is greater that 11, | | 

41 19 : BLSS LARGE _ARGD 3; branch to large argument logic | 

: ; Logic for positive medium sized arguments. Get pointer into ATAND_TABLE. | 

54 54 FC 8F 9C 57° ROTL #-4, RG, RG ; R4 = index into MTHS$$SAB_ATAN table | 

54 FFFFFFOO 8F CA 5 BICL #-256, R4 : zero high order bits of~index 

54 FC15 CF44 90 5 MOVB MTHSSAB_ATANCR4], R4 ; R4 = offset into ATAND_TABLE | 

54 FDO? CF44 DE 6 MOVAL ATAND_TABLECR4J, R4 : R4 = pointer to XHI | 

Compute 2 | 

51 84 ¢4 6 MOVL (R4)+, RI ; R1 = XHI 

520 at 9 865 MULF3 RI, Rb, R2 : RQ = XexXHl | 

2 8 40 D 66 ADDF #1, R 3 R2 = 1 + X*XHI | 

. ae) ee 69 67 SUBF Ri, RO ; RO = X = XHI 

50 52 46 b $8 DIVF R2, RO 3 RO = Z = (X = XHI)/(1 + X*XHI) | 

08 20 > Evaluate 2*P(Z*#2) | 

ZE 50 00 0406 B78 ‘ MOVL RO, -(SP) ; Push Z onto the stack | 

50 50 44 9 7 MULF RO, R : = Z**2 

Fo92 cF «02s 50s 55 0 2 POLYF RO, #ATANDLENI-1, ATANDTABI olieahs | 

: = 

50 gF 44 f 16 MULF (SP)+, Re 3; RO = ATAND(Z) = 2*Q(Z2*#2) 

50 4 40 1 7 ADDF (R4)+, R 5 Re = ATAND_XHI_LO + ATAND(Z) 

50 64 40 18 78 ADDF (R4), RO ; RO = ATANDTX) = ATAND_XHI_HI + 

| 

| 
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2-005 ATHSATAND RG - Special ATAND routine att 1 90:47 MTHRTL.SRCJMTHATAN.MAR; 1 ° Gy 
418 879 : (ATAND XHI_LO + ATAND(Z)) | 
05 0418 0 RSB ; Return 
tie BB: 
0088 31 O41C § SMALLD: BRW SMALL_ARGD ; Dummy Label used to avoid adding 
4iF 4 3; an extra insrtuction in the 
7} 5 3 medium argument logic 
41F bs? ; Large positive argument logic. 
41F H : s 
41F 8 
O41F 890 LARGE_ARGD: 
54 0000C080 8F 50 47 O41F 91 piney RO, #-1, R4 3: RG = -W = -1/XK 
4 546 45 0427 B36 MULF R4, R4, RO ; RO = Wee2 
FO73 CF 46002 S50) 55 (040B 89 POLYF RO, #ATANDLEN1-1, ATANDTABI 
0431 894 : RO = P(wee2) | 
50 54 44 0431 895 MULF R4, RO ; RO = -ATAND(Z) = <-Z*P(W#*2) 
50 00004384 8F 40 tN $38 ADDF #LF_90, RO ; RO = ATAND(X) = 90 = ATAND(Z) 
05 0438 89 RSB 3; Return 
043C 898 
043C 899; 
043C 900 ; Logic for negative arguments 
ott 901 ; | 
043C 90 
043C 903 NEG_ARGD: | 
54 50 BECO BF AS 0430 904 SUBW3 = #*XBECO, RO, R4 ; Argument is less than 3/32, | | 
19 0442 905 BLSS smal ARGD ; branch to small argument logic 
54 O36F BF Bl 0444 906 CMPW #*XO35F, RS ; Argument is greater than 11, 
~~ rh sdf BLSS N_LARGE_ARGD ; branch to large argument logic | 
Rete +t ; Logic for negative medium sized arguments. Get index into ATAND_TABLE. | 
54 54 FC 8F 9C 0448 911 ROTL #-G, R4, RG ; R4 = index into MTH$$AB_ATAN table | 
54  FFFFFFOO 8F CA 9530 aig BICL #-256, ; clear high order (unused) bits of ind 
54 FBAB8 CF44 90 045 91 MOVB MTHSS$AB_ATANCR4], R4 ; R4 = offset into ATAN_TABLE 
54 FCOA CF44 ODE Be28 aie MOVAL ATAND_TABLECR4], R4 : R4& = pointer to XHI 
046 916 ; Compute Z 
Be6 917 ; 
51 8s DO 046 918 MOVL (R4)+, RI 3 R1 = XHI 
52 50 51 45 0466 919 MULF3 R1, RO, R2 : R2 = X*XHI 
52 08 52 43 046A 920 SUBF R2, #1, R2 3; R2 = 1 = X*XH] = 1 + X*(-XHI) | 
50 351 40 pegs 921 ADDF R1, RO ; RO = X + XHI = X = (=XHI) | 
50 52 46 0471 9 ¢ DIVF R2, RO ; RO=2 
474 923; 
rhe 3 : 3 Evaluate Z*P(Z**2) 
43 50 DO 0474 9 $ MOVL RO, -(SP) ; Push Z onto the stack 
5 50 44 0477 9 MULF RO, R ; RO = Z#*2 
FD24 CF 02 50 55 Q47A 928 POL YF RO, #ATANDLENI-1, ATANDTAB1 
4 929 3; RO = P(Z**#2) | 
28 Hi 44 04 930 MULF (SP)+, RO 5 Re = ATAND(Z) = Z*P(Z**#2) 
0 4 rf : 4 931 SUBF (R4)+, RO ; RO = AT XHI_LO + ATAND(Z) 
50 64 42 04 3 9 ¢ SUBF (R4), RO ; RO = ATANDTX) = ATAND_XHI_HI + 
4 9 : (ATAND_XHI_LO + ATANB(Z)) 
05 : 4 : RSB ; Return 


EEUU EERE 


MTHSATAN 
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CJMTHATAN.MAR; 1 1) 
§ ; Logic for large negative arguments 


" 
3 N_LARGE_ARGD: 
4 DIVE 


54 00005980 5 0 47 rh I RO, #-1, RS 3; RO = W = 1/IK! 
4 $2 4 MULF R4, ; RO = Wee2 
FDO8 CF a9 946 POLYF RO, #ATANDLEN1=1, ATANDTAB1 
94 3 RO = P(Wee2) 
29 54 44 ace MULF R4, RO 3 RO = ATAND(W) = WeP(Wee2) 
50 00004384 8&F 42 45 SUBF - #LF_90, RO ; RO = ATAND(X) = ATAND(W) = 90 
0 SA 308 RSB 3; Return 
4A 94 
4A7 9948 gs 
4A7 =—-949 = Small argument Logic. 
O4A7 950 ; 
4A7— 951 
4A? 93¢ SMALL_ARGD: 
54 50 39 4A7 895 MOVF RO, R4 ; R4 = argument = X 
Hs 1 O4AA 954 BEQL $ 
50 8000 8F AA Boat 955 BICW #*x8000, RO : RO = IX! 
50 3A00 8F 861 0481 956 CMPW #*x3A00, RO 3 Compare. ‘gies to ix! 
08 19 0486 957 BLSS 1 3; Needs P 
50 54 FCFC CF 45 0488 958 MULF3  FI_OV_180_M_64, R4, RO ; RO = xeEb 4/180 - 64) 
0c 11 Q4BE 959 BRB 2$ 
0460 960 
50 50 44 pace 961 1$: MULF RO, RO 3 RO = X#"2 
FCE7 CF 02 50 55 04C 96 POL YF RO, #ATANDLEN2-1, ATANDTAB2 
0409 96 3 RO = Q(X*#2) 
50 54 44 0409 96 MULF R4, RO : RO = X#Q(X##2) 
54 0300 8F AO O4CC 965 2S: ADDW #*x300, R4 7 RG = Xe2e26 
50 54 40 04D1 966 ADDF R4, RO 3; RO = ATAND(X) = X#2e*6 + X#Q(Xe#2) 
05 0404 967 3$: RSB ; Return 
0405 968 
0405 969 
0405 970 
0405 971 END 
| 
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AIPLUS 00307 1 | 
AIPLUSD 5 R 1 
A2PLUS 11k 1 
A2PLUSD OF R 1 
CMASK = 1 
ATANDLEN1 = 0 | 
ATANDLEN2 = | 
ATANDTAB1 900 1A4 R 01 
ATANDTAB2 0000180 R 01 
ATAND_ TABLE 900000F R 01 
Aras 5 £85888 
= 
TAN Soogooe RO 
TANT 000000F0 R 01 
ATAN_TABLE 0000003C R 01 
F 00000214 R 01 
FD 000003A2 01 
LARGE _ARG 00000297 R 01 
LARGE ~ARGD 0000041F R 01 
180 = 00004434 
F-90 = 0000438 
LF~M90 = 0000¢3B4 
LF“MPI_OVER_2 = OFDBCOC9 
F7P = OFDB4149 
LETPI_OVER 2 = OFDB40C9 
LF-PI-OVER~2_HI= OFDB40C9 
LF-PI-OVER~2-LO= BD2EB4 


BD 
L = 00000004 
MTHSSAB_ATAN 00000004 
MTHS$AB_ATAN V 00000000 RG 
MTHSSJACKET_AND eeeeeree x 


MTHS$$SIGNAL geeeeene 


x 
MTHSATAN Q00001BC RG 


THSATAND 0000034A RG 
MTHSATAND2 0000035A RG 
MTHSATAND_R4 38 RG 
MTHSATAND R4D 00000 R 
MTHSATAN_R4 00000230 RG 
MTHSK_INJARGMAT teeeeeee 


pele lelelelelelelelelelelelejleajleleoleleloleola) 
Ok a ak td a oh 


~N MmMorwoe 
DDDDDDADADD ZD 
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Aw weno ren oman mamas > 


! Psect synopsis : 


poem enscea sree wose ee 


PSECT name Allocation PSECT No. Attributes | 
| 


00000000 < -} 69 ( 0.) NOPIC USR CON CL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
( 1237.) O1¢ 1. CL = SHR 


- ABS. ABS L 
-MTHSCODE 0000040 ) PIC USR CON REL L EXE RD NOWRT NOVEC LONG 


ten oece cesses sewn ne wee nea 


: Performance indicators H 


Phase Page faults CPU Time Elapsed Time 


Initialization 29 08:80 00-10 00:00:00.31 
Command processing 110 0:00:00.54 se 
Pass 109 bs FE 38 00:00:09.4 
Symbol table sort 0 00:00:00.0 00:00:06.0 
Pass 2 220 BF 28838 -00 BF OO Og 
Symbol table output § Ba: 0:0 +88 0:00:00.22 
Psect synopsis output 0:00:00.0 00:00:00.06 
Cross-reference output Sb So 88 + 8 399-09 
Assembler run totals 479 00:00:05.30 00:00:22.3 


The working set Limit was 1050 pages. : 

15739 bytes (31 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 52 non-local and 8 Local symbols. 
1031 source Lines were read in Pass 1, producing 21 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 


$ooone wom meow eo oceeooooenn + 
H Macro Library statistics H 

Macro Library name Macros defined 

-$255$DUA28: CSYSLIBISTARLET.MLB;2 0 


0 GETS were required to define 0 macros. 
There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/D1I SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHATAN/OBJ=OBJ$:MTHATAN MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) *MSRCS$: 


0257 vaxvums va ENT CORPORATION 


V4.0 3 . ND PROPRIETARY | 


